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A STUDY OF OXIDATION OF HYDROGEN BASED ON FLASHBACK OF 

HYDROGEN- OXYGEN-NITROGEN BURNER FLAMES 

By Burton D. Fine 
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SUMMARY 

The flashback of hydrogen- oxygen-nitrogen flames was studied as a 
function of pressure, burner diameter, equivalence ratio, and oxidant 
strength. The results were treated on the assumption that the product of 
the critical boundary velocity gradient for flashback and the initial 
concentration of that reactant which is not in excess is proportional to 
a mean reaction rate associated with the flame zone. It was further 
assumed that this reaction rate can be expressed in terms of initial con- 
centrations and flame temperature. 

Measurements at constant flame temperature yield orders of reaction 
with respect to hydrogen and oxygen. These do not vary with flame tem- 
perature. Measurements in which pressure is varied for several values 
of oxidant strength at constant equivalence ratio yield a total order of 
reaction and a function describing the dependence of the mean reaction 
rate on flame temperature. The total reaction order is independent of 
flame temperature and equal to the sum of the orders for hydrogen and 
oxygen. The dependence of the reaction rate on flame temperature cannot 
be described by a constant activation energy. The activation energy 
obtained apparently increases with flame temperature. Flashback results 
can be described by a single rate constant which is independent of equiva- 
lence ratio. Values were estimated for this rate constant as a function 
of flame temperature. 


INTRODUCTION 

The thermal concept of flame propagation predicts that the critical 
boundary velocity gradient for flashback should be related to the mean 
reaction time within the flame and should be independent of transport 
processes occurring through the flame zone (refs. 1 to 3). Previous 
studies have indicated that the exponent describing the dependence of the 
critical gradient on pressure is closely related to an over-all reaction 
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order within the flame (refs. 3 and 4). However, if one assumes that the 
dependence on flame temperature of the mean flame reaction rate can be 
described by an activation energy given in the Arrhenius form, the over- 
all reaction order shows considerable variation with change of diluent 
or oxidant strength for a single fuel, or on substitution among chemical- 
ly similar fuels (ref. 5). 

The present study is an attempt to evaluate, through suitable flash- 
back measurements, the various parameters in the rate expression for the 
hydrogen- oxygen-nitrogen flame as a function of flame temperature. In 
this way one might investigate, in detail, the apparent change of over- 
all reaction order with flame temperature. The hydrogen-oxygen-nitrogen 
flame was chosen for several reasons. First, the flame reaction is prob- 
ably the simplest among active flames involving fuel and oxygen. Second, 
the fact that all constituents are diatomic (and thus have similar heat 
capacities) makes possible many changes in proportion of constituents of 
the initial mixture without significant change in flame temperature. 

Third, flashback could be obtained over a range of pressures and burner 
diameters compatible with existing apparatus. 




THEORETICAL BACKGROUND 


The Semenov equation expresses the burning velocity as a function 
of rates of heat and mass transport and average rates of heat release 
occurring within a flame. The derivation of this equation and the assump- 
tions involved are presented in reference 6. According to reference 7, 
the Semenov equation may be approximated, under certain conditions, as 


tt 2 2a b% Pb 

U b T 


( 1 ) 


(All symbols are defined in the appendix.) Here ct^ is the thermal dif- 
fusivity associated with the flame zone, defined as X^/cp The 


quantity 0*13 is a mean flame reaction rate defined as 


p t 


o>. 


dT/(T, - T ) . 


The quantity a 0 may be defined by considering equation (20) of refer- 
ence 6, which relates the highest temperature attained in the flame zone 
T-^ to the amount of heat Q released by the chemical reaction. The equa- 
tion is given as 


a o Q = p o c p,o( T b - V 


( 2 ) 


Here Q is in molar units, but Cp Q is the specific heat per unit 
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mass, assumed to be constant. Thus, if T 0 is constant and T-^ is a 
maximum for nearly stoichiometric proportions of initial constituents, 
a Q is roughly equal to fuel concentration for lean flames and oxygen 

concentration for rich flames. For stoichiometric flames a D may refer 
to either fuel or oxygen, provided that Q is expressed in the proper 
units. The proportionality between a Q Q and T b - T Q further requires 

that the extent of dissociation at temperatures near T^ is not too 
large . 

The thermal quenching equation of reference 8 has the form 

_ FGa q a o( p q/ p o^ (3) 

a q — 1 ~ 

0) q 

Here the subscript q generally denotes conditions in a nearly quenched 
flame, and a D has the same meaning as before. Since o^/aiq is very 
nearly constant for similar flames (ref. 7), equations (l) and (3) may be 
combined and simplified to give 


(4) 

d q a 0 

The critical boundary velocity gradient for flashback gf- is generally 
proportional to U b /d q (refs. 2 and 3). Thus, equation (4) may be written 


¥o K “b 

One may assume that the mean reaction rate for a hydrogen- oxygen flame 
reaction may be expressed in terms of a general mass-action law 


% = k ° a H 2 a 0 2 a x F ( T b) 

where l, m, and n are constants. The quantity F(T b ) is a function of 
flame temperature alone. It is apparent that equation (6) may refer to 
the rate of disappearance of either hydrogen or oxygen or to the rate of 
formation of water. At constant flame temperature equations (5) and (6) 
may be combined to give 


m-1 n l 
gf 0(1 a H 2 a 0 2 a x 


(?) 
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for lean flames and 


g f 06 a H 2 a 02 la x (S) 

for rich flames. For stoichiometric flames, either equation vould apply. 
This result follows from the definition of a Q as the concentration in 
the initial mixture of the reactant which is not in excess. It also as- 
sumes that concentrations in the initial mixture can "be used to describe 
a reaction taking place within the flame zone at flame temperature. That 
is, it does not account for changes in concentration due to heating or 
depletion of reactants caused by reaction in the flame zone. (To avoid 
the use of double subscripts, it is noted here that all concentrations in 
equations (5) and (6) and in all subsequent discussion refer to initial 
concentrations, unless otherwise noted.) The effect of heating may be 
incorporated into FfTj,), but not the effect of reactant depletion. Ex- 
pressions for correcting initial concentrations as a function of flame 
temperature and equivalence ratio are developed in reference 6. However, 
from results of calculations of burning velocity as a function of equiva- 
lence ratio presented in reference 9, it appears that the correction sug- 
gested in reference 6 will often overcompensate for depletion. The pro- 
cedure followed here is to evaluate parameters appearing in the rate ex- 
pression in terms of initial concentrations. Before conclusions are 
drawn based on these values, however, the effect of correcting initial 
concentrations for depletion is examined. 

The form of equations (6) to (8) suggests the following methods 
based on flashback by which the various unknown parameters might be 
evaluated: 

(1) At equivalence ratios less than unity, one might replace nitrogen 
by oxygen in various proportions at a constant concentration of hydrogen. 
This substitution should maintain nearly constant flame temperature. At 

constant pressure, then, g f ex a £ , provided that 2=0, that is, that 

u 

the reaction rate does not depend on the concentration of inert gas. In 

that case, a plot of log g^> against log 3 q should have a slope equal 

2 

to n. 

(2) At equivalence ratios greater than unity, one might replace 
nitrogen by hydrogen in various proportions at a constant concentration 
of oxygen. This substitution should also maintain flame temperature 
nearly constant and, at constant pressure, should permit determination 
of the order with respect to hydrogen, again provided that 2=0. 

(3) The over-all reaction order r is defined as 


r = 2 + m + n 


( 9 ) 


) OT 
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The relation between r and g f may be shown by expressing equation (7) 
in terms of equivalence ratio as follows: 

gf oc ( 2<p) n a^ m 1 a^F(T|- ) ) (lO) 

This equation applies for lean or stoichiometric flames. If is 

varied at constant pressure and equivalence ratio , equation (10) may be 
written 


s f ^ a H 2 a x F ^ T b^ ^ 1;L ^ 

If the reaction rate depends only on hydrogen and oxygen, 1 = 0 and 
equation (ll) becomes 

g f “ a H 2 lF(T b ) (12) 

If a form for F^) is assumed, for example, the Arrhenius form with an 

assumed activation energy, equation ( 12 ) can be used to determine r in 
terms of the pressure exponent on g^ and known properties of the initial 

mixture. However, if a form cannot be assumed for F^), then equation 
(12) contains two unknowns, r and F(T^). In order to evaluate these 
two quantities, two independent flashback measurements are necessary. 

One suitable measurement would be the variation of gf with pressure 

for a single initial mixture. Another measurement would be the variation 
of gf at constant pressure and equivalence ratio with increasing pro- 
portion of oxygen in the oxidant mixture. By performing both measurements 
over a range of flame temperature, one could examine any change in r 
with T^ . Furthermore, comparison of r values so obtained with values 

obtained independently for m and n would serve as a check on the 
assumption that 1=0. 


APPARATUS AND PROCEDURE 

The low-pressure combustion chamber and the system for metering and 
mixing gases are described elsewhere (ref. 5). The main burner was a 
tube about 125 centimeters long which had a diameter of 1.89 centimeters. 
It was water-cooled near the lip. Burner diameter was varied by using 
inserts with inside diameters of 0.928, 1.016, and 1.459 centimeters. 

For flashback measurements a stable flame was established first. Then 
pressure was increased slowly at constant mass flow until the flame 
flashed back. Wherever it was possible, flashback was obtained with 
normal laminar flames. However, where the initial mixture contained 
29.6 percent hydrogen at equivalence ratios less than 1, flames were un- 
steady and polyhedral, and flashback was obtained under turbulent 
conditions . 
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Tank hydrogen (98 to 99 percent H 2 ) was used as fuel. Oxidant mix- 
tures having various nitrogen- oxygen ratios were generally prepared by 
metering the flow of tank nitrogen and oxygen separately and allowing 
the two gases to mix. Two oxidant mixtures (one 35 percent oxygen and 
65 percent nitrogen, the other 50 percent oxygen and 50 percent nitrogen) 
were specially prepared and available in tanks. All materials were used 
without further purification. 


TREATMENT OF DATA 
Reduction of Data 

For laminar flames g^ was computed from 


g f = SU f /d 


(ref. 10) . For turbulent flames the expression 

u f 0.8 

g f = 0.023 — Re f 


(13) 


(14) 


was used (ref. 10) where Reynolds number of flashback is defined as 

Re f = U f dp/p (15) 


The density was computed on the assumption that the gases were ideal. 

The viscosity was estimated from plots of viscosity as a function of 
composition (ref. ll) for hydrogen-nitrogen and hydrogen- oxygen mixtures. 


Effect of Burner Diameter 

If g f is to serve as a measure of a time associated with the flame 

reaction, it must be independent of burner diameter.. In the present ex- 
periments, this condition held for all burners used, provided the oxidant 
strength ^ was less than 0.5. For t > 0.5, however, critical flash- 
back gradients for the two smallest burners (1.016 and 0.928 cm in diam. ) 
were considerably smaller than those obtained with the two largest burn- 
ers. This effect did not depend on the material of which the burner was 
made or the mechanical condition of the burner lip. 

A review of past flashback experiments performed with the present 
experimental arrangement shows that this deviation follows a trend pre- 
viously observed. Thus, for flames having burning velocities in the 
range 200 to 300 centimeters per second (hydrogen-air and propane- 
oxygen-nitrogen mixtures with ^ = 0.5) burners 1.89 to 0.55 centimeter 
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in diameter gave consistent flashback gradients; a burner 0.3 centimeter 
in diameter caused flashback at much lower values of For flames 

having burning velocities of about 400 centimeters per second (propane- 
oxygen and ethylene-oxygen-nitrogen mixtures with ^ = 0.5), a burner 
0.55 centimeter in diameter tended to deviate to lower values of gf* 
Thus, the present deviation by a 1-centimeter burner seems to be as- 
sociated with burning velocities of about 600 centimeters per second or 
larger. 

The cause of this decrease in g f is not presently known. The 
deviation is favored by high burning velocities and small diameters; it 
is probably peculiar to the present apparatus, since such an effect has 
not been previously reported. The procedure followed here is to ignore 
those results obtained in 1-centimeter burners which show deviation 
toward small gf values. 

Other slight deviations, associated with low Reynolds numbers, are 
shown as dotted lines in figure l(a). These had been observed previous- 
ly (ref. 3) and attributed to partial quenching of the flame by the 
cooled wall. 


RESULTS AND DISCUSSION 


General Description of Results 


Critical boundary velocity gradients for flashback are shown as 
functions of pressure within the combustion chamber in figures 1 and 2. 
Results are shown for various values of burner diameter and initial com- 
position. Since the initial gas mixtures contain three components in 
generally varying proportions, two quantities must be specified to de- 
fine the composition. Those chosen are the equivalence ratio 


<P 



stoich 


and the oxidant strength 


(16) 


= 



( 17 ) 


Values of for stoichiometric hydrogen-air flames shown as dashed 

lines in figures l(c), l(d), and 2 are based on previously reported 
measurements , which were obtained near an equivalence ratio of 1 
(ref. 4). 
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Thus, the dashed lines represent an interpolation between cp = 0.95 and 
<p = 1.10 for laminar flames and (p = 0.95 and (p = 1.20 for turbulent 
flames (ref. 12). 

All flame temperatures discussed are calculated equilibrium tempera- 
tures obtained by the method of reference 13. They are based on a single 
set of thermochemical data and are therefore consistent among themselves. 
In general, slopes with respect to pressure have been estimated by eye. 

If least-square methods were used, certain points would have had to be 
arbitrarily weighted and others disregarded, particularly at low Reynolds 
numbers, where quenching may have affected the data, and at high Reynolds 
numbers, where instability in the flow may have affected the results. 
Accordingly, the second decimal place in the pressure slopes probably has 
little physical significance. However, it was retained through subsequent 
calculations . 


Results at Constant Flame Temperature 

Results near 2580° K. - As discussed in the section THEORETICAL 
BACKGROUND, the individual contributions of fuel and oxygen to the over- 
all order of reaction can be determined from the change in the critical 
flashback gradient with change in the fuel or oxygen concentration in 
the original mixture at constant pressure and flame temperature. The ap- 
proach outlined in that section has been applied in the following way. 

A stoichiometric flame having an oxidant strength of 0.35 and a calcu- 
lated flame temperature of about 2560° K at a pressure of 15 centimeters 
of mercury was chosen as an experimentally convenient reference. The 
change of g f with pressure for this flame was determined. Then, 
similar measurements were made on mixtures obtained from the reference 
mixture by replacing nitrogen by oxygen at constant hydrogen concentra- 
tion. All such lean mixtures have flame temperatures equal to about 
2560° K at a pressure of 15 centimeters of mercury. One may also replace 
nitrogen by hydrogen at constant oxygen concentration. If the reference 
stoichiometric flame is again taken to be that for which \|r = 0.35, all 
rich mixtures obtained by such substitution will have flame temperatures 
calculated to be about 2620° K at 15 centimeters of mercury. Thus, 
values for the order with respect to fuel and oxygen determined in this 
way refer to flame temperatures which differ by about 60° K. 

Flashback results for flames at these temperatures are shown in 
figures l(a) and (b). The variation in flame temperature with pressure 
over the range of pressures for which data were obtained is less than 
the difference between calculated temperatures on the lean and rich sides 
and need not be considered. On the lean side (fig. 1(a)) the concentra- 
tion of hydrogen is kept constant at 41.0 percent. On the rich side 
(fig. 1(b)), an oxygen concentration of 20.5 percent is maintained. Since 
d In g f /d In P is constant, within experimental uncertainty, for 
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lean, stoichiometric, and rich flames , the temperature difference of 
60° K between the lean and rich side is not significant. The fact that 
d In gf/d In P is constant permits cross plots of against fuel con- 

centration for lean flames at constant pressure with reasonable assurance 
that results so obtained will be independent of pressure. Such cross 
plots are shown on the right in figures l(a) and (b). It is seen that 


log gp is linear in log or log X- 


V 


The fact that (5 log gf/d log 



is constant indicates that 


the dependence of the reaction rate on the concentration of inert gas is 
zero. By logarithmic differentiation of equation (7) at constant hydro- 
gen concentration, one obtains 


. n + LM (is) 

P log a 0 2 / a H 2 v lo e a 0 2 / ^2 


The value of (S log a x /d log aQ )a^ 

2 2 


is sufficiently large and suf- 


ficiently variable that 7 must be very close to zero for the exponential 
dependence of gj. on oxygen concentration to appear constant. A similar 

argument, based on equation (8), applies for the apparent constancy of 
m. The results of the cross plots at constant pressure give n = 0.92 at 
pressures of 15 and 11 centimeters of mercury and m = 1.37 at 13 centi- 
meters of mercury. The calculations are summarized in table I. The sum 
of m and n is about 2.29, which agrees rather well with values for 
the over-all reaction order for the hydrogen-air flame obtained from 
measurements of flashback, burning velocity, and quenching distance 
( refs . 4 and 14 ) . 


Results near 2300° K. - It is of interest to repeat the measurements 
described in the preceding section under conditions yielding significantly 
different values of flame temperature. This should give some information 
regarding the behavior of individual and over— all reaction orders with 
changes in flame temperature. It was not feasible to work at flame tem- 
peratures much above 2600° K; therefore, measurements were made at lower 
flame temperatures. The stoichiometric hydrogen-air flame, having a cal- 
culated temperature of 2348° K at a pressure of 1 atmosphere, was chosen 
as a reference. Results are shown in figures l(c) and (d) and summarized 
in table I . 


For rich flames flashback followed a normal course. A cross plot of 
log gp against log ag^ 14.8 percent oxygen and P = 19 centimeters 

of mercury gave m = 1.39, a result in excellent agreement with that ob- 
tained near 2580° K. On the lean side, however, flames of 29.6 percent 
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hydrogen were polyhedral in shape and unsteady and displayed a regular 
spinning motion. The relation between gj and P observed for such 

flames was complicated and clearly unrelated to the flame reaction. 
Therefore, the attempt was made to obtain values of m and n on the 
lean side for mixtures containing 29.6 percent hydrogen by observing the 
flashback of turbulent flames. Previous studies on the flashback of 
turbulent flames had indicated that critical boundary velocity gradients 
for laminar and turbulent flashback differ by a factor that is independent 
of pressure and nearly the same for all systems studied (refs. 3 and 12). 
Thus, one should obtain the same values of m, n, and l regardless of 
whether measurements are based on laminar or turbulent flashback. Fur- 
thermore, since turbulent flashback seems to be controlled by processes 
taking place in a nearly laminar portion of the turbulent boundary layer, 
results obtained from turbulent flashback should apply at the temperature 
of the corresponding laminar flame, even though the space-average tem- 
perature of the turbulent flame might be somewhat lower. 

Accordingly, results are shown in figure l(d) for air (interpolated 
for cp = 1 from the data of ref. 12) and for mixtures with i|r = 0.35 
and 0.50. Calculated flame temperatures for this series of lean flames 
were nearly constant, but the value for lean flames was about 60° K lower 
than that for the reference stoichiometric flame (see table I). However, 
as shown on the right side of figure l(d), all three points fall on the 
same straight line. The slope has a value of about 0.92, in agreement 
with the value for n found at the higher flame temperature. 


Comparison of results based on initial and effective concentra- 
tions. - According to the Semenov theory, the concentrations which affect 
the flame reaction rate are effective concentrations in the flame zone. 
Expressions for these effective concentrations are given in reference 9. 
However, in figure 1, flashback gradients are plotted as functions of 
initial concentrations. It is important, then, to see how the results 
would look in terms of concentrations corrected by the expressions given 
in reference 9. For a hydrogen- oxygen- inert gas flame the relation 
between a^ and a^ can be expressed approximately as 
u 2 2 



<p<l. An equation symmetrical in 


would apply for 


<p >1. The quantity 5 is a constant whose order of magnitude is given 
by (RT^/E act )/(T- b - T Q ). According to reference 6, the value of this quan- 


tity should not exceed 0.1 for the assumptions in the theory to be justi- 
fied. Therefore, it will be assumed that 5 = 0.1, even though the value 
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Of E act which follows from this assumption, about 50 kilocalories, is 
rather large compared with other values of E ac t estimated from flame 
properties . If one defines an effective mole fraction as 

X* = a* (RT b /P), one may use equation (l) to relate initial and effective 
mole fractions at constant pressure and flame temperature. For lean 
flames, one obtains 



0 . 45 Xrr 
h 2 


( 20 ) 


and for rich flames 

x h 2 = x h 2 " 1,80 x o 2 


( 21 ) 


One may evaluate Xq and 


y* 


from data given in table I. 


* * 

Figure 3 shows plots of log against log Xq^ and log 

under the same conditions as in figures l(a) and (b). Results in terms 
of X* are considerably more complex than in terms of X. Botham ana 
n are considerably less than 1 and are strongly dependent on X . Con- 
sequently, no simple relation exists among m, n, and r. There is no. 
prior reason for rejecting results in terms of X* in favor of those in 
terms of X. However, if one accepts results in terms of X as vali , 
one cannot use them to develop a conventional mass-action expression. 

In order to develop such an expression from flashback data, one must 
assume that the expressions for effective concentrations given in refer- 
ence 9 overcorrect- for depletion of reactants and that initial concen 
trations are a more satisfactory approximation to true effective 
concentrations . 


Results at Constant Equivalence Ratio 


Figure 2 shows g f plotted against chamber pressure at a constant 
equivalence ratio of 1 for varying values of oxidant strength.. The . 
slopes of the lines increase slightly but significantly with increasing 
ilr. This increase is of the same order of magnitude as that previously 
observed for hydrocarbon- oxygen- inert gas flames (ref. 3). One may re- 
late the slopes to the reaction order r by differentiating Ration 
(12) logarithmically with respect to pressure at constant <t> and 


d In g f 
d In P 


= r 


1 + 


a m F(T b ) d m T b 
a in Tb d In P 


( 22 ) 
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As a first approximation on F(T b ) one may assume that the depend- 
ence of the rate of the flame reaction on flame temperature is a simple 
exponential function which can be described by a constant activation 
energy. For this model, which is discussed in references 3, 4, and 14, 
one may write equation (22) in the form 


d ln Sf n E act ( 1 _ , 1 \ d ln T b 

d In P " r 2R \T b + T q J d ln P 


(23) 


and evaluate an approximate reaction order r'. For consistency with 
previous applications of equation (23), an activation energy of 23 kilo- 
calories per mole was assumed for the hydrogen- oxygen flame reaction. 

The flame temperature was related to T , the temperature of a nearly 

quenched flame, by 


T q = 0.8 T b + 0.2 T q (24) 

(ref. 7). Values of d ln T b /d ln P are based on flame temperatures 
calculated at pressures of 1 and 0.1 atmosphere. 


Values of r' based on equation (23) are shown in table II. It 
can be seen that the apparent reaction order increases with flame tem- 
perature and that the value of r' for t = 0-35 is significantly 
greater than the value of the over-all order previously determined as 
the sum of orders for hydrogen and oxygen. These two factors suggest 
that the assumption of a simple exponential dependence of the flame re- 
action rate on flame temperature Is inadequate, at least for the hydrogen- 
oxygen flame reaction. 


However, if F(T b ) is treated as an unknown function, the data 
plotted in figure 2 can be used along with values of the approximate 
reaction order r' to determine the form of F(T b ) and evaluate its 

dependence on flame temperature. With this information, equation (22) 
can be used to compute a better approximation to the over-all reaction 
order at various flame temperatures. The procedure is outlined in the 
section THEORETICAL BACKGROUND and is based on the fact that, as shown 

by equation (12), gp/ajj is a function of flame temperature alone if 
equivalence ratio is held constant. A plot of log (gp/^gg ) against 


l/T b is shown in figure 4 for a constant pressure of 15 centimeters of 
mercury and a constant equivalence ratio of 1. This quantity should be 
proportional to F(T b ). The plot shows considerable curvature, which 


supports the previous suggestion that the dependence of the reaction 
rate on flame temperature could not be described by a constant activa- 
tion energy, at least for the present flame system. One may obtain 
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approximate values of d In F(T b )/d( l/T^) for various values of T b from 

the slopes of tangents drawn to the dotted curve in figure 4. The actual 
slopes must be multiplied by -0.434/T-^ to give 5 In F(T^)/S In T-^ . 
These values may be used in equation (22) to recompute a set of better 
approximations to r. The results of this calculation are shown in the 
first column on the right in table II. The values obtained scatter be- 
tween 2.15 and 2.34; they show no significant trend with flame tempera- 
ture, in contrast with the set of r' values. Of course, the accuracy 
of this calculation is not very high, since it depends on evaluations of 
derivatives from a curve which is itself not especially well defined by 
experimental data. However, the results suggest that, within the rela- 
tively large experimental error, the reaction order has a value of about 
2. 2 or 2.3 independent of flame temperature. This agrees well with the 
value of about 2.3 obtained from the determination of individual orders 
at constant flame temperature. One may use this value in equation (12) 
and compute a second approximation to F(Tp ) ). This is shown as the solid 
line in figure 4 . The two sets of data points define curves which are 
similar in shape and position. Specifically, the curvature in the plot 
of log F(T b ) against l/T^ is not sensitive to assumptions made re- 
garding the reaction order. 

It should be noted that results for hydrogen-argon-"air" flames 
(ref. 3) shown in table II and plotted in figure 4 are in satisfactory 
agreement with values for corresponding hydrogen-oxygen-nitrogen flames. 
This is a further indication that the nature of the diluent does not 
affect the course of the hydrogen oxidation at these temperatures. 


Comparison of Results for Lean and Rich Flashback 
at Constant Flame Temperature 

For the several curves plotted in figures l(a) and (b), both oxygen 
concentration and equivalence ratio are varied. Therefore, the data 
plotted in these figures may be further cross-plotted as a function of 
equivalence ratio at constant pressure. Such a cross plot is shown in 
figure 5 for a pressure of 13 centimeters of mercury. The curve consists 
of a lean and a rich branch and has a minimum at an equivalence ratio of 
1. The lean branch corresponds to a flame temperature of about 2540° K, 
whereas the rich branch represents a temperature of about 2620° K. The 
curve as a whole is not significantly different from that which would be 
obtained for a single flame temperature of about 2580° K. Thus, it is 
analogous to plots of mass burning velocity against equivalence ratio at 
constant flame temperature that are presented in references 15 and 16. 

The fact that the curve in figure 5 goes through a minimum at an equiv- 
alence ratio of 1 reflects the fact that the concentration product portion 
of the rate expression in terms of initial concentrations goes through a 
minimum for a stoichiometric initial mixture. 
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Figure 5 shows that g^ is a double-valued function of equivalence 
ratio. If g£ measures a reaction time, then one should he ahle to 
describe lean and rich branches by the constants x^ and x r , respec- 
tively, which are related to the specific rate constant for the reaction 
between hydrogen and oxygen. On the lean side one may assume a relation 
of the form 


s f\ - x i a S 2 a S 2 


(25) 


and on the rich side 




f 0. 


“ X rV a O, 


m n 
2 u 2 


(26) 


The constants x^ and x r are most easily computed from equations (25) 
and (26) by use of the cross plots shown on the right in figures 1(a) and 
(b). Thus, in terms of mole fractions equation (25) may be solved for 
x^ to give 


S* 


h - 


xS' 1 x; 




r(S 

0 O 


m+n-1 


(27) 


and equation (26) becomes 


= 


«f / goV 

c m x n_1 \ p / 


m+n-1 


X„ XL 

H 2 °: 


(28) 


Using any convenient sets of g_ values from figures 1(a) and (b) , one 

a I C ..... 


obtains 
great as x 


= 9. 1*10 C 


Thus x, is about twice as 


and x = 4. 3*10 . 

r l 

. From the data of figures l(a) and (b) it is easily seen 
that the product a^j for a given value of g f on the lean side 

equals a2 a n for the same value of g on the rich side at constant 

h 2 o 2 i - 

pressure and flame temperature. Thus, the observation that x^/x^ is 

about 2 is consistent with the fact that the steady rate of consumption 
of hydrogen, which is proportional to gf a H 2 .> i s twice the rate of 

consumption of oxygen, which is proportional to Sp a o • 


M 
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One may define a rate constant "by the expression 

k(T) = x/3 


(29) 


where 3 is a constant. The variation of k with T should he repre- 
sented by the best approximation to F(T^) given in figure 4. It is 
apparent that in order to evaluate k(T) one must know 3 or must be 
able to calculate 3 from a known value of -da^/dt, the average rate 

of consumption of hydrogen in the flame. A calculation of -daj^/dt is 
made in reference 7 for a stoichiometric hydrogen-air flame at a pressure 
of 1 atmosphere. This calculation, based on burning velocity and quench- 
ing distance, gave a value of 169 moles per liter per second. The value 
of gx. for a stoichiometric hydrogen-air flame is known to be about 
8500 second'*- 1 - at 1 atmosphere (refs. 4 and 10) . Therefore, one may use 
the relation 


r d \ 

Sf^g - " P dt 


(30) 


to evaluate 3* A value of 0. 60 is obtained. Since 3 is a constant 
independent of temperature and pressure, one may use it in equation (29) 
along with the value for determined in the previous section to ob- 
tain the result k(2580) = 1.46*10^ (moles-*-* ^) (liters -*-‘^)(sec ). If 

the dependence of the reaction rate on flame temperature is given by the 
best approximation to the quantity proportional to F(T^) shown in fig- 
ure 4, then one should be able to write 


k(T) = k°F(T b ) (31) 

If one solves equation (31) for k° in terms of k(2580) and F(2580) 
one obtains k® = 1* 55*10^. One may then plot k(T) as a function of 
T^. Such a plot is shown as a solid line in figure 6. 


It is possible to obtain values of k over a range of flame tem- 
perature by a different method. First, one may evaluate k for a 
stoichiometric hydrogen-air flame by 


- da H 2 

dt 


1.38 0.92 


= k(2348)a H Sq 


(32) 


This gives k(2348) = 0. 66xl0 7 . One may then extrapolate each of the 
curves in figure 2 to a pressure of 1 atmosphere and estimate values 
of -dap^/dt by assuming a relation of the form 
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'sM atm) a " 

. H2 . 

2 

g f (l atm)a H ^’ 

1 


(da /dt) 

_2 2,1 atm 

(da^ /dt) 

n 2 1, 1 atm 


(53) 


where the subscripts denote any two values of oxidant strength. Inser- 
tion of resulting values of -da^/dt into equation (32) gives k as a 

function of T^. These results are shown as the dotted curve in figure 6. 


The agreement between the two curves is satisfactory, considering 
the many possibilities for systematic error involved in their determina- 
tion, It is difficult to decide which of the two is preferable. The 
dotted curve involves several long extrapolations of low-pressure data 
up to 1 atmosphere. The solid curve assumes that (3 is independent of 
flame temperature, an assumption which cannot be checked by use of 
existing data. 




SUMMARY OF RESULTS 

The flashback of hydrogen-oxygen-nitrogen flames has been measured 
at reduced pressures as a function of burner diameter, pressure, equiva- 
lence ratio, and oxidant strength. On the basis of a treatment which 
assumed the critical boundary velocity gradient for flashback to be a 
measure of a mean flame reaction time, the following conclusions were 
drawn: 

1. The order of the flame reaction with respect to hydrogen is about 
1.4, independent of flame temperature. 

2. The order of the flame reaction with respect to oxygen is about 
0.9, independent of flame temperature. 

3. The total order of reaction is independent of flame temperature 
and equal to the sum of the orders with respect to hydrogen and oxygen. 

4. The variation of flame reaction rate with flame temperature 
cannot be described by a constant activation energy. 

5. From relations among flashback, burning velocity, and quenching 
distance, a specific rate constant for the flame reaction can be esti- 
mated as a function of flame temperature. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, October 1, 1958 
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APPENDIX - SYMBOLS 
concentration, moles/liter 
specific heat, cal/(g)(°K) 
burner diameter, quenching distance, cm 
activation energy, kcal/mole 

dimensionless coefficient, (Tq - T Q )/(T b - T 0 ) 

dimensionless geometric factor 

critical boundary velocity gradient, sec 

specific rate constant, (moles )(liters )(sec ) 

exponential dependence of reaction rate on concentration of inert 
constituents, dimensionless 

exponential dependence of reaction rate on hydrogen concentration, 
dimens ionle s s 

exponential dependence of reaction rate on oxygen concentration, 
dimensionless 

pressure, atm or cm Hg 

heat of reaction, cal/mole 

gas constant, cal/ (mole) (°K) or (liters) (atm) /°K 
Reynolds number, dimensionless 
total reaction order, dimensionless 
temperature, °K 
time, sec 

mean stream velocity, cm/sec 
burning velocity, cm/sec 
mole fraction 

thermal diffusivity, cm /sec 


dimensionless coefficient 
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6 dimensionless additive term 

x rate constant in terms of flashback 

X thermal conductivity, cal/ ( sec) (cm) (°K) 

[i viscosity, poises 

p density, g/cm^ 

(p equivalence ratio, dimensionless 

\|/- oxidant strength, dimensionless 

cn reaction rate, moles/ (liter) (sec) 

Subscripts : 
b normal flame 

f flashback 

H 2 hydrogen 

Ng nitrogen 

0 2 oxygen 

o initial conditions 

q quenched flame 

r,Z rich and lean equivalence ratios 

x chemically inert constituents of original mixture 

Superscripts: 

— mean value 

0 independent of temperature 

1 first approximation 

* corrected for depletion of reactants 
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TABLE II. - DETERMINATION OF OVER-ALL REACTION ORDER 
FROM FLASHBACK OF STOICHIOMETRIC HYDROGEN- 


OXYGEN- NITROGEN FLAMES 


- — 

Oxidant 

Mole 

Flame 

d In 

d In g f 

Total 

strength, 

+ 

fraction 

of 

temperature, 

d In P 

d In P 

reaction 

order 


hydrogen, 

Xtt 

H 2 

% 

1 atm 



r ! r 

0.21 

j 0.296 

2348 

0.017 

1. 33 

2.25 2.26 

a . 21 

a . 296 

a 2590 

a . 027 

a l. 51 

a 2. 38 a 2. 31 

.28 

.365 

2570 

.025 I 

1. 54 

2.41 2.34 

.35 

.411 

2683 

.030 : 

i 1.55 

2.41 2.17 

. 50 

. 500 

2834 

.034 | 

; 1. 72 

2. 57 2. 15 

. 70 

.583 

2936 

.040 ; 

1. 79 

2. 62 2. 16 


Nitrogen in air replaced by argon or helium. 
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Figure 1. - Flashback of hydrogen-oxygen-nitrogen flames at constant flame temperature. 





Figure 1. - Continued. Flashback of hydrogen-oxygen-nitrogen flames at constant flame 
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Figure 1. - Continued. Flashback of hydrogen-oxygen-nitrogen flames at 
constant flame temperature. 






Figure 1. - Concluded. Flashback of hydrogen-oxygen-nitrogen flames at constant flame 
temperature . 
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Figure 2. - Flashback of stoichiometric hydrogen-oxygen-nitrogen flames 
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First approximation 






Figure 5. - Dependence of critical boundary velocity gradient for flashback on equivalence ratio at constant temperature and 
pressure. Flame temperature, 2580° K; pressure, 13 centimeters of mercury. 









i ’“I c * 

d ^ £ , 
3 • 3 

*3 -p ‘ 


CO 


CO 


co 


U 

co 


•o = 

II 

mS 


g | | g w 

O 5 O 3 .S < 



< 

z 


co ° 
QJ c 
p p 


•'-' 0) 

« u _ 

O O ’ 

g § ^ 

> 0 ) 5 

hS 2' 


qj 


QJ . 


^ O 

o * 

0) d 


3 u 
QJ 3 
bfl CO 
>» co 
X o> 

0 «-• 

1 a a; 

QJ 13 3 
bD ® 01 
O g CO 
t-, p QJ 
T3 T3 U 
>, a) a, 
-a ^ 


o O T3 
O C 
_ JO QJ 
O * P« 

ia« 

C« > QJ 

g ? J3 
*-» c H 

j- QJ 

5 bo 

> £ QJ 

-US 

^ d 


CU QJ 


CO QJ 

0> & 

3 .. 

rt g 
> p 

-p 2 

« 3 

2 O 

I ° 
, 0} 

o g 


« co 

^ nt 


d qj 

£ -P 


o * 5 : 

^ "2 QJ 

2 a 


o fl 
a) d qj~ 

5 w S 

-t_» d 

cl ^5 Cl 

XJ O 
-*-* d d 

§S.s 

•2 co £ 


d P +j 
u .2 d 


d -n 0> 

g e 6 

p d H 


3 


d 

d 'O 




S H 2 ■S 

qj 3 •*•; 

a c g S 

g o -a -a i 

d rj m a> i 

a rrt Cl Q<j 

2 ■» « < 

"O ,, 


-M t-« g 
3 u .2 
d a +j 

111 

° g 8 

QJ H Y 

*3 « 

Sh *h 


d 


co 


q> T3 


:13- 


Z Z < & 


H ** jo i 


„ c e 

2 g bo . 
* 2 >» rr 
> * * 


.5 P 


qj d 

If 

-M QJ 

si 

<u h 


3 
Qj , 


W 

co 

LfJ 


s 

d 1 

J 


i .o LT 5 


CO 

co 

l 3 2, _ co 

h a ^ 


CO 


5 a 
n s 

3 g 
CQ S 


.o . - < 

E 

o 3 .S < 

U In I* Z 


»-H csj »_ J_3 


|| 

*S ! 

;-i 

3 | 

z ! 

E ! 

£ i 


QJ 


QJ > -O 



d ! 


- p o 

p 2 5 

5 -S 

w 03 S 

g £ QJ CO 

o 3 .2 < 
O Om Pm Z 



co ° 

QJ C 

S 3 

3 o 
c 

C 3 
QJ ^ 


H s 

> 0) w 

>> S 2 


■8 


i QJ 

C Sh 


QJ ' 

£3 5j * 
u g co 

c ^ c 
O - 
5 bo 


1 T 3 
) ci « 


QJ 


tH co < 


3 *rt 

bo co 2 

>. co d 

X QJ 

0 ^ - 

1 a QJ 
C _. U 
QJ 13 3 
bo ® cn 

O y CO 
i~> P QJ 
T 3 ^ t-< 

>. g a 

•P 

O •* s 
'■! « » 


5 

n 

.2 71 


o ^ 


Jh 'S 

3-3 


« ■> 

2 

g- * ^ 

SI V 


2 a 

o a - 

Q) d qj 

5 w a 

4-* ‘ d 
rt HirJ 

o 

*- rt 

c 2 


3 tn 


o 

p 

2 * 
o -*-* 

d -rj 

“ C 


* lH 


CO ^ 
1 _ , 


d 

JO 


CO * ’P 

a a - 


d 

2 ® 

o 

3 

c 5 

.2 * 


rt 3 2 

QJ Qj r- * 

X 


rj 

a d 

QJ 


c5i.2 


£ < 


H ^ JO 


3 ^ u 
O bo d 


dj qj 

C bC 

s >» 

* x 

Z o 


a 3 ■“ 

2 d o 

to > 3 

^2 

H .2 P . 

,l s g 

AJ M 

g ^ u « 

g 

S - S e 
S g *5 

f » a g 

1 ga| 

CD ° g rt 

3 QJ g ^ 

d * P 

^ 53 * 

2 o 3 g 

O S ° § 

c u _: ° 

^ a| 

3 « rt ^ 
QJ qj 3 

a^ p -o 

QJ t H o 
■ to 

o ^ 

QJ g U 
u p co 
3 g QJ 
QJ ^ T 3 


3 


’g m 

p o 


d « 3 * 2 , 




fill 

O U 


Q - 

3 O 

2 g 
3 W 
CQS 

3 < 

Z. 

~ a 


< 

z 



o 



